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We present our views on the issues raised in the chapter by Griffin and Zaremba [A. Griffin and E. Zaremba, 
in Quantum Gases: Finite Temperature and Non-Equilibrium Dynamics, N. P. Proukakis, S. A. Gardiner, M. J. 
Davis, and M. H. Szymanska, eds., Imperial College Press, London (in press)]. We review some of the strengths 
and limitations of the Bose symmetry-breaking assumption, and explain how such an approach precludes the 
description of many important phenomena in degenerate Bose gases. We discuss the theoretical justification for 
the classical-field (c-field) methods, their relation to other non-perturbative methods for similar systems, and 
their utility in the description of beyond-mean-field physics. Although it is true that present implementations of 
c-field methods cannot accurately describe certain collective oscillations of the partially condensed Bose gas, 
there is no fundamental reason why these methods cannot be extended to treat such scenarios. By contrast, 
many regimes of non-equilibrium dynamics that can be described with c-field methods are beyond the reach of 
generalised mean-field kinetic approaches based on symmetry-breaking, such as the ZNG formalism. 



I. INTRODUCTION 

In Ref. Hi, Griffin and Zaremba (GZ) have offered a cri- 
tique of the so-called classical-field (c-field) method, compar- 
ing it to the Beliaev broken-symmetry approach to Bose su- 
perfluidity [2 3 1, and formalisms derived on this basis, such 
as the ZNG method HHJ. A large part of their discussion 
concerns the validity of the separation of the Bose field oper- 
ator into coherent (classical) and incoherent parts in practical 
applications of the classical-field approach, which they con- 
trast with the separation into condensed and non-condensed 
parts in symmetry-breaking theories. Regarding the applica- 
tion of these methods to the description of experimental sys- 
tems, GZ focus on the collective modes of oscillation of finite- 
temperature BEC, which the ZNG formalism has reproduced 
with great success, and reason that the classical-field approach 
is unsuitable for the description of such phenomena. 

Here we explain why the absence of a spontaneous 
symmetry-breaking assumption is not a failing of the 
classical-field approach, but is actually an important feature 
which ensures its generality and utility in describing many 
interesting scenarios of degenerate Bose-gas dynamics, in- 
cluding non-quasistatic growth of the condensate [6], regimes 
of matter- wave turbulence Q, and low -dimensional sys- 
tems (8]|9]. We also argue that the success of the ZNG method 
in modelling the collective oscillations of finite-temperature 
condensates does not appear to be due to many of the formal 
consequences of the symmetry-breaking assumption cited in 
Ref. (T]. Although present implementations of the classical- 
field method lack the dynamical description of the above- 
cutoff atoms required to accurately reproduce some collective 
oscillations, there does not seem to be any fundamental reason 
why it cannot be extended to include these dynamics. How- 
ever, in our view the classical-field method is very useful even 
without such an extension, as it is the only existing method- 
ology which can treat the non-equilibrium dynamics of the 
system in strongly fluctuating regimes. 



H. SYMMETRY BREAKING 

A major point identified in Ref. [lj is the absence of an 
a priori distinction between condensed and non-condensed 
modes, such as that which accompanies symmetry breaking, 
in the classical-field approach. In our view, this is not a fail- 
ing of the method. As GZ note, the defining characteristic of 
Bose superfluidity is the emergence of a new thermodynamic 
variable on the superfluid side of the phase transition. The 
new variable that arises is, fundamentally, the superfluid ve- 
locity, which underpins Landau's famous two-fluid model of 
superfluidity IflOl . This velocity field is in general attributed 
to the presence of an underlying condensate [11 1, and is de- 
termined by the gradient of the condensate phase. A defi- 
nite value for the phase [U(l) symmetry breaking] is thus not 
required to define the superfluid velocity; this was noted by 
Anderson Ifl2l . who nevertheless advocated a description of 
superfluidity based on symmetry breaking, as it provides a 
unified description of both superfluidity and macroscopic in- 
terference phenomena in condensed systems [ 13 1. As is well 
known, true spontaneous symmetry breaking only occurs in 
infinite systems [14], and for an isolated system of atoms, the 
particle-number superselection rule formally prohibits the ap- 
pearance of a finite field expectation value 11131 . Never- 
theless, field-theoretical calculations are typically formulated 
in a 'restricted ensemble' in which a finite first moment of the 
field does exist [16]; such a restricted ensemble can be intro- 
duced by adding symmetry-breaking terms to the Bose-field 
Hamiltonian 0, which correspond to an external phase ref- 
erence and thus lift the superselection rule (see Ref. ifTTl and 
references therein). In this chapter we do not discuss the fun- 
damental validity of this approach, but concern ourselves only 
with the pertinent operational qualities of symmetry-breaking 
theories. It is of course true that many profound and ele- 
gant consequences follow immediately from the assumption 
of Bose symmetry breaking 0. GZ note in particular that the 
symmetry-breaking approach makes it evident that the veloc- 
ity field of the condensate exactly determines the local veloc- 
ity of the total superfluid density — which is in general not 
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equal to the condensate density — and yields a rigorous def- 
inition of this superfiuid density lfl8l . However, it does not 
appear to be necessary to appeal to symmetry breaking to de- 
fine these quantities in general 1 19 1. 

Moreover, the natural relationship between the condensate 
and the superfiuid density in the symmetry-breaking approach 
is of little consequence in the application of many formalisms 
based on symmetry breaking, such as the ZNG method J4] [5) . 
GZ stress that the symmetry-breaking separation of the sys- 
tem into condensed and non-condensed (thermal cloud) parts 
in the ZNG method gives a 'natural way of capturing the 
two-fluid nature of superfiuids resulting from an underlying 
Bose condensate'. However, in existing implementations of 
the ZNG method, the non-condensate atoms are treated in 
a single-particle (Hartree-Fock mean-field) description, and 
as such, the superfiuid component is precisely the condensed 
component of the system (described by a generalised Gross- 
Pitaevskii equation). One would expect that the implemen- 
tation of a Bogoliubov-quasiparticle model of the thermal 
cloud (Chapter 7 of Ref. [5|) would describe the small non- 
condensed contribution to the superfiuid density in the com- 
paratively straightforward case of a three-dimensional con- 
densate. However, scenarios of experimental interest in which 
there is a pronounced difference between the condensate and 
the superfiuid density, such as the two-dimensional Bose gas 
— in which thermal fluctuations erode long-range order but 
leave the superfiuid density relatively unaffected 120 -221 — 
cannot be treated in the ZNG approach, precisely because 
such systems do not conform to the limiting case of a distinct 
condensate with well-defined quasiparticle excitations. 

It is important to note that there are many such features of 
the condensed Bose gas which are fundamentally beyond a 
description in terms of a symmetry-breaking approach. The 
symmetry-breaking ansatz assumes that the amplitude of the 
condensate orbital is a classical variable that, as GZ note, does 
not undergo any fluctuations, and cannot exhibit any higher- 
order correlations with the non-condensed component of the 
field. By contrast, as one approaches the transition to the 
normal state the condensate population exhibits increasingly 
large fluctuations ll23Tl26l . which require a description using 
non-perturbative methods [27 1, and more generally particle- 
number conservation implies anticorrelation between the con- 
densed and non-condensed populations E51 . 

GZ also review the consequences of the assumption of sym- 
metry breaking for the excitations of the system. In par- 
ticular, they refer to the fact that the Bose broken symme- 
try leads to the condensate and non-condensate components 
sharing the same single-particle excitation spectra, and that 
the spectrum of density fluctuations of the system is also 
locked to the single-particle spectrum BJ. These are in- 
deed some of the most profound and elegant results of the 
symmetry-breaking approach. However, approximate field 
theoretical formalisms based on symmetry breaking almost 
uniformly fail to reproduce these formal results. This is the 
famous Hohenberg-Martin dilemma |[T6l l29l : all approxima- 
tions for the self energies in the conventional Green's func- 
tion approach either violate conservation laws (and thus do 
not yield an acoustic density-fluctuation spectrum), or vio- 



late the Hugenholtz-Pines theorem lfl6l |30"1 . which requires 
a gapless single-particle spectrum. As a result, the equiva- 
lence of the single-particle and density-fluctuation spectra is 
never realised in these approaches, and this is also true of the 
ZNG formalism. By contrast, the so-called dielectric formal- 
ism QTl . an alternative perturbative approach which evades 
the Hohenberg-Martin dilemma and ensures coincidence of 
the two spectra, has not produced an accurate model of high- 
temperature collective oscillations [32|. The undeniable suc- 
cess of the ZNG method in treating collective oscillations of 
finite-temperature condensates |5, 33-35 1 therefore does not 
appear to be related to many of the formal features of the 
symmetry-breaking field-theoretic framework cited by GZ. 

Finally, we note that although it does seem necessary to in- 
voke symmetry breaking to provide a basis for a general field- 
theoretic approach to perturbation theory in the presence of 
a condensate [16, 36], the predictions of symmetry-breaking 
methods for the excitation spectrum, including the Bogoli- 
ubov, Hartree-Fock-Bogoliubov, and Beliaev-Popov approx- 
imations, are also obtained in explicitly number-conserving 
approaches (see references J37] ED, ED, and EQl 52, re- 
spectively). 



HI. BASIS OF THE CLASSICAL-FIELD 
APPROXIMATION 

The classical-field method is distinguished by the sepa- 
ration of the Bose field operator into two components: a 
part composed of single-particle modes which comprise the 
coherent region (or condensate band), and the remaining 
modes that constitute the complementary incoherent region 
(see Ref. [42]). GZ contrast this distinction between the 
two regions with the division of the field operator into con- 
densed (classical mean field) and non-condensed parts in the 
Beliaev method. They note that the classical field v Fc( r ) m 
the classical-field approach has (in contrast to the symmetry- 
breaking condensate of the Beliaev approach) 'no special sta- 
tus' and that the classical-field separation is 'artificial'. It is 
indeed true that the separation of the Bose field into coherent 
and incoherent parts is not based on a strict physical distinc- 
tion, and as such, v Fc( r ) is not to t> e interpreted as being in 
any fundamental way physically distinct from the remainder 
of the Bose field as, for example, the condensed component 
of the field is in the Beliaev approach. Nevertheless, this sep- 
aration is an important component of the c-field method. 

The division of the field in the c-field method can be mo- 
tivated as follows. The lowest-energy components of the in- 
teracting Bose field cannot be understood in terms of single- 
particle-like modes or excitations, as interactions strongly 
couple the single-particle modes of the corresponding non- 
interacting system (e.g., the harmonic-trap eigenstates). In the 
limit of a well-defined condensate, this strong coupling is pre- 
dominantly representative of the fact that the excitations of the 
condensed gas are Bogoliubov quasiparticles, rather than, e.g., 
dressed single-particle (Hartree-Fock) states 11431 . In more 
general situations, such a quasiparticle picture may be inap- 
plicable, and the coupling of modes may be better understood 
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in terms of density and phase fluctuations of an equilibrium 
quasicondensate (in low dimensional systems) [44-48], or as 
a regime of strongly turbulent behaviour, such as is predicted 
to occur in the late stages of non-adiabatic condensation ||49l . 
Moreover, near the superfluid phase transition the field ex- 
hibits critical fluctuations |50|, which herald the breakdown 
of perturbation theories and preclude the interpretation of the 
system in terms of a well-defined condensate and quasipar- 
ticle excitations in this regime BUI IBTI . In all of these sce- 
narios, the low-energy portion of the field is characterised by 
correlations beyond a simple Hartree-Fock model of dressed 
single-particle states. 

Fortunately, it is also the case, at least for moderately high- 
temperature regimes of degenerate Bose gases, that the single- 
particle modes that span the part of the field that exhibits non- 
trivial correlations are classically occupied; i.e., they have 
mean occupations {N^} » 1, such that their quantum fluc- 
tuations can reasonably be neglected. As such, although the 
correlations and dynamics of this region may be quite non- 
trivial, the effects of quantisation of the low-energy modes are 
comparatively subdued; i.e., these correlations and dynamics 
result primarily from the multimode, self-interacting character 
of the field, not from its nature as a quantum field. 

At higher energies, the occupations of field modes subside, 
so the highest-energy modes cannot be treated in a classical- 
field approach, and the excitations of the system also revert to 
a single-particle-like structure [52 1. Provided that the return of 
the excitations to an essentially single-particle nature occurs 
at a lower energy scale than that at which the classicality of 
mode occupations becomes violated (i.e., at which (N^) < 1), 
one can introduce a division of the system into two parts (see 
Ref. l53l ). such that one part contains all the non-trivial cor- 
relations and dynamics, whereas the other contains all modes 
which have sub-classical occupations. One thus identifies a 
component *i*c(r) of the system which is potentially beyond 
any mean-field treatment, but which is, however, amenable to 
a classical-field description, while the complementary high- 
energy component 'i'i(r) can be treated in a simple mean-field 
approach 11541 . 

The separation of the two components in the c-field ap- 
proach is therefore indeed an artificial one, designed to al- 
low for the application of the classical-field approximation to 
modes for which it is required, while not spuriously applying 
it to weakly-occupied modes for which it is not valid. This 
discussion gives us guidelines as to the limits of the cutoff 
in the projected Gross-Pitaevskii equation (PGPE) / stochas- 
tic PGPE (SPGPE) approach: the cutoff should be at least of 
order /j ~ g\(p\ 2 above the condensate eigenvalue, in order to 
include all quasiparticle structure and critical fluctuations in 
the classical field, while being not greater than k B T above the 
condensate eigenvalue in order to ensure reasonably classical 
mode occupancies ((Nt) > 1) [42| (see also Ref. Il55l0 . Within 
these limits the choice of cutoff is arbitrary, and results should 
be practically independent of the precise choice of cutoff in 
this range l54l . 

It is important to note that there are two distinct, but 
closely related, classes of classical-field models. Microcanon- 
ical classical-field methods such as the PGPE involve closed- 



system, Hamiltonian equations of motion for the classical re- 
gion of the Bose field, and neglect any coupling to the com- 
plementary incoherent component. By contrast, the grand- 
canonical SPGPE [ 56 1 formalism — which unifies the PGPE 
approach with the quantum kinetic theory of Refs. [52, 571l58ll 
— and the closely related Stoof stochastic Gross-Pitaevskii 
equation (SGPE) [59| involve explicitly stochastic equations 
for the classical field, with noise and damping terms which 
represent the effects of coupling to the above-cutoff region. 
Although implementations of these methods to date assume 
the above-cutoff region to be in thermal equilibrium even 
when studying system dynamics Il6l l60l466ll . this is not a fun- 
damental restriction of such formalisms |59, 671 1681 . 



IV. EQUILIBRIUM CORRELATIONS IN THE 
CLASSICAL-FIELD APPROACH 

As the classical-field approach is not based on a division 
of the field into condensed and non-condensed components, 
quantities such as the condensate must be inferred from the 
field correlations a posteriori [8 69 1. In cases where a well- 
defined condensate exists, it can be found by applying the 
Penrose-Onsager definition of condensation [70| to the clas- 
sical analogue of the one-body density matrix [69|. More gen- 
erally, properties of the system such as quasicondensation and 
superfluidity can also be deduced from correlation functions 
of the field. 



A. Microcanonical Classical Fields 

The extraction of equal-time correlation functions from 
the PGPE is straightforward at equilibrium, as the nonlin- 
ear Hamiltonian evolution of the c-field samples the corre- 
sponding microcanonical density defined by the conserved 
energy and other first integrals of the system 11421 . Corre- 
lation functions of the field, such as the one-body density 
matrix, are thus formally given by averages over this mi- 
crocanonical density, which can be approximated by aver- 
ages over time. It is important to note that, having extracted 
the Penrose-Onsager condensate orbital from a classical-field 
simulation, one can define a 'fluctuation field', composed 
of the part of the field orthogonal to the condensate orbital, 
but taking into account the relative phase between this field 
component and the condensate ll26l . This field is then the 
c-field analogue of the number-conserving non-condensate 
field operators introduced in symmetry-preserving quantum- 
field approaches 13714401 I7T1 1721 . In this way, one can cal- 
culate [in addition to the coherent-region part of the non- 
condensate density «'(r)] the so-called anomalous moments 
of the non-condensed component of the field [26] in the 
U(l)-symmetric microcanonical ensemble of the PGPE (see 
also Refs. 117314751 ). These correlations correspond to many- 
body processes, which yield corrections to the condensate- 
condensate and condensate-non-condensate scattering ll40l 
17614781 and, in particular, to the chemical potential of the con- 
densate. An analysis of a three-dimensional condensate at 
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finite-temperature equilibrium [26] shows that the anomalous 
averages can be significant, and moreover that the conden- 
sate obtained from a Penrose-Onsager analysis of the classi- 
cal field appears to be consistent with the mean-field picture of 
the condensate as a nonlinear eigenfunction of a generalised 
Gross-Pi taevskii operator that involves pair and triplet anoma- 
lous moments ||4l[76) (see Ref. [74] for a related study of the 
pair anomalous average with the one-dimensional SGPE). 

The extension of the ZNG method to describe the non- 
condensed component of the gas in terms of Bogoliubov 
quasiparticles (see Chapter 7 of Ref. |5|) would — as in the 
generalised mean-field kinetic treatments of Refs. 17914811 — 
include the pair anomalous average, and the attendant cor- 
rection to the condensate's self interaction. We note, how- 
ever, that in present implementations of the ZNG method 
the anomalous average is neglected, despite not necessarily 
being small a priori. Moreover, it is not so clear how the 
higher-order correlations of the non-condensate fluctuations, 
such as the triplet correlations, could be built into the ZNG 
model. The presence of these correlations in c-field equilib- 
ria underlines not only the naturally self-consistent nature of 
the Penrose-Onsager condensate in the classical-field equi- 
librium, but also the fact that the c-field method furnishes a 
non-perturbative description of the field fluctuations. 

We also note that as the PGPE is a fundamentally dynamic 
equation of motion, it (as well as other Hamiltonian classical- 
field methods [82 83 1) also gives access to (approximate) dy- 
namic correlations at equilibrium [73 83-90], as in the 'Lan- 
dau dynamics' approach to spin models 19"T1 . Importantly, 
this dynamical character is not merely introduced artificially 
to provide a means of sampling the thermodynamic ensem- 
ble (cf. Ref. [92 1), but forms an approximation to the actual 
dynamics of the Bose field: for large mode occupations, the 
evolution of the bosonic field is well described by a classical 
field equation, from which the quantum-field dynamics can be 
approximated by higher-order corrections in the inverse mode 
occupation [93 1. 

B. Grand-Canonical Classical Fields 

The stochastic GPEs l56l l59l describe the c-field as a 
stochastic process, with correlations obtained in principle 
from averages over an ensemble of distinct trajectories. In 
regimes where the underlying Hamiltonian part of the evolu- 
tion behaves ergodically, one expects that a single trajectory 
will cover the appropriate thermodynamic ensemble densely, 
as in the usual Langevin equation approach to sampling ther- 
mal distributions [50]. One can therefore substitute time av- 
erages for ensemble averages at equilibrium. However, in 
low-dimensional regimes, one must be mindful of the possi- 
ble proximity of the underlying Hamiltonian part of the evo- 
lution to nearby integrable models (see, e.g., Ref. l94l ) as, 
in general, the addition of noise and damping terms may not 
overcome the tendency of the system to remain 'trapped' in 
quasi-regular regions of phase space [95 1, preventing efficient 
exploration of the thermodynamic ensemble. Correspond- 
ingly, studies of one-dimensional equilibrium systems with 



the Stoof SGPE explicitly consider an ensemble of distinct 
trajectories El ED E2 EH ED . 



C. Fluctuations 

An important feature of the classical-field approach, noted 
by GZ, is its description of field fluctuations. In general, it 
is only in the limit of a well-defined condensate that a clear 
division into a mean-field condensate and a complementary 
non-condensate part is valid. As noted in Section[ll] the over- 
all amplitude of an otherwise well-defined condensate may 
undergo fluctuations, which become important as the system 
approaches the transition to the normal phase. More gener- 
ally, and particularly in low-dimensional systems, the concept 
of a condensate may not be useful in describing the system, 
and the appropriate 'mean-field' approaches in such cases take 
the alternative route of describing the system in terms of den- 
sity and phase fluctuations about an equilibrium quasiconden- 
sate Il44l447ll . As the classical-field method is not based on 
the assumption of the existence of a well-defined condensate 
or quasicondensate, it is equally applicable to those regimes. 
In particular, direct comparisons between the Stoof SGPE and 
the modified Popov theory of Andersen et al. l46l WT\ have 
shown good agreement between the two (for related studies 
see Refs. 1191 1741196^991 ). In general, one expects that a true 
condensate, if present, is given by the part of the field which is 
both density- and phase-fluctuation suppressed l47l l48l, and 
SGPE equilibria have been found to be consistent with this 
identification ll74l[T00l . 

An important related question raised by GZ is to what ex- 
tent a classical-field approach can provide an adequate de- 
scription of the critical region near the phase transition. It is 
well known that the transition to Bose condensation belongs 
to the universality class of the classical XY model If23"lll011 . 
As such, any classical |</<i 4 -type model, such as that consid- 
ered in the PGPE, describes the universal characteristics of the 
Bose-condensation transition, and such classical models are 
commonly used to investigate the critical physics of the Bose 
gas 11101141031 . As such approaches treat the system without 
recourse to any assumption of spontaneous symmetry break- 
ing, the issues of infrared divergences alluded to by GZ (see, 
e.g., Ref. ED ) do not arise (though such divergences are nat- 
urally cut off in the experimentally relevant case of harmonic 
trapping 11041 ). The new feature of the PGPE classical-field 
approach as applied to the harmonically trapped case, is that 
one expects, having properly taken account of the trapping po- 
tential, to also include non-universal features peculiar to the 
system being studied. Provided that the cutoff is chosen ap- 
propriately, the method should yield quantitative predictions 
for the system's critical behaviour, and indeed, the predictions 
of the PGPE approach for the T c shift of the harmonically 
trapped, interacting Bose gas agree with experiment to within 
the present experimental uncertainty l54l . 

In light of this discussion, it is potentially confusing that the 
classical component of the Bose field operator in classical- 
field methods is sometimes described as an 'order parame- 
ter' ll59l l68l . Since by construction *Fc(r) describes a range 
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of low-energy modes in addition to the condensate (should 
one be present), we agree that it is incorrect to think of 
^cW as a formal order parameter (except when the cut- 
off is low enough for this to coincide with the condensate 
1 105]). However, appears as a field described by an 

effective time-dependent Ginzburg-Landau equation (partic- 
ularly in the stochastic GPEs |56, 59 1), and it is unfortunately 
conventional in the relevant literature to refer to the field ap- 
pearing in such descriptions as an order parameter [ 106 1; this 
should of course not be confused with the order parameter as- 
sociated with the Bose-condensation transition. 



D. Superfluidity 

GZ note that the symmetry-breaking approach relates the 
superfiuid flow directly to the velocity field of the underlying 
condensate, and that this feature does not have to be inserted 
into the symmetry-breaking theory as an additional assump- 
tion. We stress, however, that such an assumption is not re- 
quired in implementing the classical-field method. In general, 
the superfiuid density is given by the response of the field to 
an applied velocity field (or phase twist), which leads to var- 
ious correlation functions from which the superfiuid density 
can be calculated 1107141 10l . The superfiuid density can 
be extracted from the c-field using these theoretical expres- 
sions l98l , just as (for example) the superfiuid density of the 
XY model is calculated from classical Monte Carlo calcula- 
tions II 1 II . 

It is true that each of the eigenmodes of the one-body den- 
sity matrix calculated in the classical-field approach defines its 
own individual velocity field; but this would of course be true 
of any one-body density matrix of a Bose system, obtained in 
any approximation. As GZ note, one expects that inasmuch 
as the system exhibits a well-defined condensate, it is this par- 
ticular mode which defines the superfiuid velocity field lfl9l . 
From the velocity fields of all the eigenmodes of the one-body 
density matrix, one can construct a 'hydrodynamic' velocity 
field, which corresponds to the total mass current |112|. This 
hydrodynamic velocity field bears, of course, no a priori re- 
lation to superfluidity. However, in a system which contains 
a superfiuid component which is significantly larger than the 
condensate — i.e., one in which there is a significant contri- 
bution to the superfiuid density from non-condensed atoms — 
the associated supercurrent should constitute a contribution to 
the hydrodynamic current that matches the velocity field de- 
fined by the condensate orbital, and preliminary c-field simu- 
lations of quasi-two-dimensional superflows suggest that this 
is indeed the case II 131. 



V. NON-EQUILIBRIUM DYNAMICS IN THE 
CLASSICAL-FIELD APPROACH 

A. Non-Equilibrium Correlations 

An important focus of Ref. [1] is the viability of the 
classical-field method as a means of describing the non- 



equilibrium dynamics of Bose superfluids. GZ correctly note 
that although the substitution of time averages for ensemble 
averages is reasonable for equilibrium systems, it is not ap- 
plicable in non-equilibrium situations, in which one must in 
general explicitly consider averages over an ensemble of dif- 
ferent classical-field trajectories. These are distinguished by 
distinct choices of initial conditions for the field (and in the 
stochastic GPEs by distinct samples of the dynamical noise 
processes |114|), which may be sampled from an equilib- 
rium thermal distribution [6, 115|, or from the Wigner dis- 
tribution corresponding to the Bogoliubov vacuum in simula- 
tions starting from zero temperature ll85ll . Although averag- 
ing over an ensemble of trajectories may in general lead to 
situations where there is no well-defined condensate, it does 
nevertheless seem at least reasonable to interpret individual 
classical-field trajectories as representative of individual ex- 
perimental realisations of the system ll42l l64l [TOBI . The indi- 
vidual trajectories, while having no direct formal meaning in, 
e.g., the truncated Wigner interpretation, certainly closely re- 
semble the images seen in experiments. Moreover, aside from 
constructing /orma/ correlation functions, the ensemble of tra- 
jectories does allow one to amass summary statistics for, e.g., 
decay rates of solitons [64 66] and vortices l65l . and frequen- 
cies of topological-defect occurrence ]6] II 16L By contrast, 
mean-field approaches such as the ZNG method only describe 
a 'mean' trajectory for the condensate in non-equilibrium evo- 
lution 111711118 11. and give no information on the possible sta- 
tistical dispersion of such quantities in experiments. 

In quasi-equilibrium dynamical regimes, one may find that 
the timescale on which the macroscopic dynamics of the con- 
densate take place is sufficiently long, compared to the co- 
herence timescale of thermal fluctuations, that one can un- 
ambiguously extract a non-equilibrium condensate and ther- 
mal cloud from a single classical-field trajectory by averaging 
over suitable intermediate timescales Il88llll9l . However, no 
such obvious separation of timescales appears in general non- 
equilibrium scenarios, making it difficult, if not impossible, 
to identify a single, well-defined condensate mode during the 
evolution. We emphasise that this is, in some sense, precisely 
the point: in general non-equilibrium situations one may not 
be justified in assuming a priori that a well-defined conden- 
sate exists. 



B. Collective Oscillations 

A particularly important regime of non-equilibrium dynam- 
ics discussed by GZ is that of collective oscillations of finite- 
temperature condensates, the accurate description of which 
has presented a profound challenge for theories of Bose- 
condensate dynamics since their first realisations in experi- 
ments 1 120] 11211 . The results of extensive theoretical stud- 
ies by many authors 11341 l68l II 22141 26l reveal that the cou- 
pling between the oscillations of the condensate and those of 
the thermal cloud is of crucial importance for describing the 
anomalous frequency up-shift of the m = quadrupole oscil- 
lation at high temperatures. A recent study of the collective 
oscillations in a classical-field (PGPE) approach [115] found 



6 



increasing down-shifts of the m — 2 oscillations with increas- 
ing temperature, consistent with previous studies 1122L but 
did not observe the high-temperature up-shift of the m — 
mode (see also Ref. 11271 ). As noted by the authors of ref- 
erence II 151 . and by GZ, this is perhaps not surprising, as 
the dynamics of the above-cutoff incoherent region were ne- 
glected in that study. 

VI. DISCUSSION 

We emphasise here that all dynamical processes within the 
coherent region are described (in the classical-field limit) by 
the c-field methods. It is only processes which couple across 
the artificial boundary to the above-cutoff incoherent region 
which are neglected in the PGPE approach. In the stochastic - 
GPE approaches, these processes are described by the dissi- 
pative and stochastic terms that represent the coupling of the 
c-field to the above-cutoff region. However, these methods 
do neglect the dynamic evolution of the above-cutoff region, 
and the influence of the c-field on these dynamics. We stress 
that making a classical-field approximation for the low-energy 
modes of the gas does not in itself necessarily preclude con- 
structing a model that includes these processes. The coupling 
between the two components is in principle provided for in 
Stoof's path-integral derivation of the SGPE U28I . and 
although the Gardiner SPGPE is derived on the basis of an ex- 
plicit 'tracing out' of the above-cutoff region ll56ll . its deriva- 
tion could be extended to include a kinetic description of the 
above-cutoff atoms. In such extensions, the above-cutoff re- 
gion could be treated in a Hartree-Fock approach (as in cur- 
rent ZNG implementations), and its contribution to the non- 
condensate mean-field potential 2gn'(r) could easily be in- 
cluded in the c-field equation of motion. (This contribution 
is typically neglected on the basis that it is flat over the co- 
herent region at high temperatures [54|, but would, as GZ 
note, be essential for the accuracy of the coupled dynamics.) 
Such an approach should ultimately yield a method which 
includes both the 'egalitarian' treatment of highly occupied 
modes in the classical region, while also describing the cou- 
pling of the superfluid (quasi-)condensate to the complete re- 
mainder of the gas; i.e., it would include all the physics of 
the ZNG method, while additionally allowing for non-trivial 
fluctuations and correlations in the low-energy component. 

However, GZ assert that the 'natural' way to include 
the coupling between the condensate and the thermal 
cloud of non-condensed atoms is by means of a many- 
body perturbation-theory approach based on the symmetry- 
breaking concept. Although it is clear that the ZNG method 
introduces the crucial features necessary to describe the par- 
ticular non-equilibrium dynamics of collective oscillations, 
we do not agree that it is therefore preferable to pursue such 



a broken-symmetry approach to the non-equilibrium problem 
in general. It is conceivable that the ZNG method could form 
the basis for a generalised method in which fluctuations of, 
e.g., the phase of the condensate about the mean-field solu- 
tion could be calculated from the fluctuation-dissipation the- 
orem (as in Stoof's approach |59|), and used to augment the 
ZNG predictions ex post facto, thereby extending the ZNG 
approach to describe more general states of the Bose field in 
(quasi-)equilibrium regimes. However, this approach would 
still be based on the a priori assumption that an underlying 
well-defined condensate exists, which may not be justified in 
more general non-equilibrium scenarios. Indeed, there are sit- 
uations in which any approach based on symmetry breaking 
is incapable of yielding a description of the field dynamics 
which is even qualitatively correct (see, e.g., Ref. 11291 ). 

In the view of GZ, the utility of the classical-field method is 
that it allows one to 'address some non-trivial non-equilibrium 
problems without having to introduce a lot of formal machin- 
ery typical of field theoretic calculations based on the Beliaev 
formalism'. We hope that this discussion has clarified why 
we do not regard the classical-field method as a mere means of 
simplifying field theoretical calculations, but as a distinct tool, 
with a broad range of applicability that is largely complemen- 
tary to that of kinetic methods based on symmetry breaking, 
such as the ZNG formalism. 
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